Towards measuring Entanglement Entropies in Many Body Systems. 



o 
o 

(N 

X) 
D 

oo 
(N 



I 

o 
o 



> 

o 
o 
cn 
o 
\o 
o 



X3 
O 

o 



X 



Israel Klich/ Gil Refael,^ and Alessandro Silva^ 

^ Department of Physics, California Institute of Technology, MC 114-36 Pasadena, CA 91125 
^ Abdus Salam ICTP, Strada Costiera 11, 34100 Trieste, Italy 

We explore the relation between entanglement entropy of quantum many body systems and the 
distribution of corresponding, properly selected, observables. Such a relation is necessary to actually 
measure the entanglement entropy. We show that in general, the Shannon entropy of the probability 
distribution of certain symmetry observables gives a lower bound to the entropy. In some cases this 
bound is saturated and directly gives the entropy. We also show other cases in which the probability 
distribution contains enough information to extract the entropy: we show how this is done in several 
examples including BEC wave functions, the Dicke model, XY spin chain and chains with strong 
randomness. 



Entanglement entropy was first considered as a source 
of quantum corrections to the entropy of a black hole 
[1] . Since then its study expanded considerably in extent 
and purpose. Today its evaluation allows an estimate 
of the effectiveness of physical systems as quantum com- 
puters, since the entanglement entropy characterizes the 
resources in our disposal to perform computations. As 
importantly, the entanglement entropy is emerging as a 
tool in the field of many body systems. Close to quantum 
phase transitions and at critical points it was shown to 
exhibit universal properties [2, 3]. In particular, in one 
dimensional conformally invariant systems the entangle- 
ment entropy is a measure of the central charge of the 
underlying conformal field theory [4]. 

Consider a system partitioned into two parts A and B. 
The entanglement entropy Se in the quantum state p = 
evaluates how many qubits in A are determined 
by (or entangled with) qubits in B. It is defined as the 
von Neumann entropy Se = — Tr[p^ \og2{pA)] of the 
reduced density matrix pA = Tre [p] ■ 

As it stands, the notion of entanglement entropy is 
rather abstract. Indeed, to the best of our knowledge, de- 
spite the large body of literature dedicated to the study of 
entanglement entropies in many body systems, the actual 
physical possibility of measuring this quantity, and so in- 
vestigating experimentally its properties (such as scaling 
behavior), has not been addressed. This is the main mo- 
tivation of this Letter. 

It is clear that, given the entire distribution of corre- 
lation functions and observables one may always recover 
the density matrix, and hence Se- A step in this direc- 
tion was taken in Ref. [5] , where it was emphasized that a 
quantum state may be reconstructed from the set of sus- 
ceptibilities with respect to a large enough class of exter- 
nal potentials via the Hohenberg-Kohn theorem. In Ref. 
[5] , this was used to calculate the linear entanglement en- 
tropy and negativity of at most 4 dimensional subspace 
(for example two spins) with a large system. It is, how- 
ever, unclear how this can be implemented to find the 
entanglement between two large subsystems where the 
entire density matrix cannot be realistically measured. 

In this letter we show that, given a quantum many 



body system, one may always identify a class of nat- 
ural observables, whose fiuctuations can be related to 
the entanglement entropy. After some general consid- 
erations, we provide several examples of systems where 
the entanglement entropy can be estimated, or directly 
extracted, from the probability distribution P{x) of the 
possible outcomes x of such observables. This connection 
will take different forms in different problems. 

A key ingredient of our analysis consists in consider- 
ing directly, instead of Se-, the so called "measurement 
entropy" S[0] of properly chosen observables O, i.e. the 
Shannon entropy S[0] = —J2x-Pi^)^'^sP{^) associated 
with the probability distribution P{x) of the outcomes x 
of O [6] . In classical systems this quantity, measurable by 
definition, is always lower then the overall entropy. On 
the other hand, in quantum systems it can in general be 
either larger or smaller than the entanglement entropy 
[7] . Below we select a class of local observables such that 
S[0] provides a lower bound on the entropy Se- While 
in some cases the bound is saturated and S[0] = Se, 
in other cases, where estimates based on S[0] do not 
capture the scaling behavior, we are nevertheless able to 
extract Se from the distribution P{x)- 

The class of observables whose measurement enables 
us to estimate the entanglement entropy is constructed 
as follows. Given a state ip of interest, let us denote by £ 
the set of observables O — Oa®I+I®Ob, acting locally 
on A and B, for which ijj is an eigenstate. Note that C is 
nonempty [8]. The Schmidt decomposition of ip can be 
written as V' = X^cf |q!,j > ®\s — a,i >, where s is the 
eigenvalue of O acting on i/', and such that OA\ct,i >~ 
a\a,i > and Ob\s — a,i >= (s — a)\s — a,i > (here 
i ranges over the degeneracy of eigenstates of Oa with 
value a) 

One may write the reduced density matrix as pA — 
irsp = '^PaPcx, where we have defined pa = 



>< a,i\ with Pa = SilcfP- Therefore, 



for the entanglement entropy one has 

Se = SpA] - P^^^Po' log P°' ^ ^[^^ 



(1) 



where S'[Oa] is the measurement entropy associated to 
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the probabihty distribution Pa. 

The inequahty in Eq. (1) is completely general. Inter- 
estingly, the equality Se = S[Oa] is realized if and only 
if tipa log Pa = for all a, as for example in the case 
of no degeneracy of the eigenvalue a, or in more inter- 
esting cases where the pa's describe pure states. Such 
systems do indeed exist: we will consider, for example, a 
BEC like state, a two mode squeezed state, and the Dicke 
model. As remarked above the measurement entropy has 
the advantage of being directly measurable. By repeated 
measurements of Oa we can recover the distribution of 
outcomes, and so extract S'fO^]. Moreover, we note that 
Se = MaxoecS[OA]- To see this, choose the operator 
Oa diagonal with non-degenerate eigenvalues Aj in the 
basis \iA) appearing in the Schmidt decomposition of ijj, 
and Ob diagonal with eigenvalues — A^ in the basis lis), 
ip is an eigenstate of eigenvalue of Oa I + I (S> Ob , 
and SfO^] = Se- Note that in the same way, if we define 
'S'l[0] = ^ P{x){l — P{x)) to be the linear measurement 
entropy of O, it obeys Sl > Sl[0], where Sl = 1 — trp^ is 
the linear entropy of the system. In the case S[0] = Se 
one also immediately has Sl[0] = Sl- 

Intuitively, a natural subset of £ is the set of "con- 
served" operators, i.e., sums of local operators which 
commute with the Hamiltonian of the system and thus 
are in C. For instance, the total spin operator for spin 
chains with rotational symmetry. These are the cases 
we consider. We remark, however, that since not for 
every Hamiltonian the set C necessarily contains such 
conserved observables [9] , in the general case, the appro- 
priate choice of O may requires a more elaborate analysis. 

Let us now illustrate the general considerations above 
with two examples in which the equality S[0] = Se is ac- 
tually satisfied. As a first illustrative example, consider a 
BEC condensate of bosons who share a particular Gross- 
Pitaevskii wave function f{x), constrained to occupy a 
box of volume V. One can imagine dividing the box into 
two parts A and B, of volumes Va and Vb respectively. 
The condensate wave function can then be written as 

\/W. 

N 

fe=0 

where N is the number of bosons, is the bino- 



Pk log2(Pfe),withPfe = Cf(«2)^-fc(l-«2)fc^ Con- 
sequently the measurement of entropy of particle num- 
ber in A in this case is exactly the entanglement entropy. 
Assuming for simplicity Vb » Va, and taking f{x) as 
uniform [f{x) = 1/^Va + Vb] then in the limit TV » 1, 
one obtains Se — 1/2 log[A''u^(l — u^)], which, in the 
thermodynamic limit (TV +oo, keeping u = N/V con- 
stant), simplifies to Se — 1/2 log(VAi^)- 

A similar example of entangled state, but this time of 
two distinct degrees of freedom, is a two-mode squeezed 
state 



4> = exp 



10 > 



= —ul E e'"* [tanh(r?)]" | | 71)2, (3) 
cosh 77 

^ " n=0 

where £, = rj e'* is the squeezing parameter, and ai, 02 
are the annihilation operators relative to the two modes. 
Identifying A and B with the two modes, tp is an eigen- 
state of the sum of local operators O = ni(E)I+I(E){—n2)- 
Since the eigenstates of ni and n2 are non degenerate 
the entanglement entropy is Se = S[hi]. In particular, 
Se = log2[l+'T-]+log2 [(1 + n)/n\n. where n = [sinh(r7)]^ 
is the average number of photons per mode. In the limit 
rj +00, one obtains S12 ~ log2[ra]. 

Squeezed states can be generated in a number of phys- 
ical situations. A particularly interesting realization is 
obtained when a collection of two level atoms (subsystem 
A) interacts with a single mode of the EM field (subsys- 
tem B): the Dicke model. In turn, the Dicke model can 
be realized in a number of ways, both using traditional 
cavity QED, as well as employing solid state circuits [cir- 
cuit QED][10]. Assuming for simphcity the photon mode 
to be in resonance with the atoms, the Hamiltonian is 



A 



H = LUoJ^+uJoa^a+-=={a^ +a) {J+ + J-), (4) 



/2j 



where J^. = ^^jlj^ [k = ±,z] are collective oper- 
ators describing the dynamics of the collection of 2j 
two level atoms. Working in the subspace where, at 

A = 0, all atoms are in the ground state, in the large 
j limit, one can conveniently use the Holstein-Primakoff 
representation Jz — b'^b — j, J+ = 6^ ^■\/2j — fe+fe^ , 

J- = (^\/2j — 6'i'6^ b, where b are bosonic modes [11]. 
In the j — > -|-oo limit, one obtains H = ujob'^b + ojoa^a -|- 



mial coefficient, and a'^,b'^ create a particle in A,B X{a'^ + a) {b^ +b) - juo- This Hamiltonian describes two 
respectively, i.e., = ^ J^^^ f {x)'ip^ {x)dx and b^ = 
7T=f/.eB/(^)V^^(^)da:,"with = l^^\f{x)\Mx. 
Note that ip is an eigenstate of the particle number oper- 
ator N = /^g^ dx'iljl.'ijjx+ J,j,^g dxijjlijjx- We may consider 
just the subspace consisting of applications of and b^ 
to the vacuum, and ignore other modes. In this subspace 
there is no degeneracy since every occupation number 
appears only once and so we have simply iS^ = 5'[a^a] = 



coupled harmonic oscillators. In particular, as a function 
of coupling A, the Hamiltonian has a Quantum Criti- 
cal Point at Ac = wo/2, where the ground state symme- 
try with respect to parity is spontaneously broken. Set- 
ting A < Ac, it is convenient to introduce the operators 

^ :(« + «^)' y = + ^^)' = - «) 

— b). Writing down the Hamilto- 



x = 



and Py = 

nian in terms of these operators, and making the rotation 
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qi = (x + y)/\/2, q2 = (.t — ?y)/\/2 one obtains a quadratic 
Hamiltonian, the spectrum being that of two independent 
harmonic oscillators of frequencies K e^^'', where K = 
^u;^ - AX-^Lul and e^" = ^(co-^ - 2Aa;o)/(a;^ + 2Aa;o). 
The ground state can be written as 

Introducing $„(a;), the eigenstates of the operator 
H^{K) = {pI + K a;2)/2, it is holds that 

-. +00 

^ = ^Mi) ^ [tanh(r?)]" $„(a;)$„(j/). (6) 

'' n=0 

Therefore, one immediately finds that, defining the oper- 
ator d = Hx{K)®I + I® {-Hy{K)), the ground state of 
the model is an cigenstate with zero eigenvalue. In addi- 
tion, Se = S[Hx{Ky\. It is interesting to notice that, for 
A <C Ac, one can approximate — loq {a)a + while 
for A ~ Ac one has ~ In analogy with the previ- 

ously considered two mode squeezed state, for A ~ Ac the 
entanglement entropy is given by <Sb ~ 2 log2(sinh(?7)) Ri 
1/4 log[p^] [11]. 

The lower bound on Se provided by the measurement 
entropy has to be considered with care, in particular in 
extracting the scaling of Se with subsystem size La in 
one [or higher] dimensional systems. In these physical 
situations, it is convenient to choose the operators O 
as corresponding to extensive observables, such as the 
total magnetization, or particle number. In this case, 
one expects the distribution of measurements outcomes 
X of O to be Gaussian in the limit of large subsystem, 
i.e. P{x) -^=i=e-(^-'^)'/2''' for La oo. For such 

observables we have S[d] = 1 + ilog(27rV< Ax^ >). 
More generally, given the variance v<Aa?>, the for- 
mula above gives always the maximal measurement en- 
tropy of a continuous variable, as one can easily checked 
by a variational argument. In other words, the variance 
provides an upper bound to the measurement entropy. 

While for such extensive operators one immediately 
obtains S[0] w Iog(\/ < Ax"^ >) < Se- in many cases 
this inequality is of limited use, since, while S[0] scales at 
most logarithmically with the variance, the entanglement 
entropy scales in fact as the variance itself and not its 
logarithm. This implies that, in many cases, the variance 
is a useful entanglement measure. Indeed the variance 
is decreased in any local physical POVM measurements 
respecting super-selection rules [12]. 

To demonstrate this, let us now consider the XY spin 
chain [4] and cases of non- interacting fermions [13]. Un- 
der the Jordan- Wigner mapping one can treat the XY 
spin chain as a special case of non-interacting fermions 
on a ID tight binding model, with no external poten- 
tials. The role of the total S% observable 5^,^^ S'f is 
played by the fermion number operator X^jg^ctcj. The 



entanglement of free fermions is particularly interesting 
as it exhibits scaling violating area law in higher dimen- 
sions [14, 15]. Indeed, a formula for the entropy scaling 
was presented in [14] in d dimensions, and was recently 
checked numerically for certain models in 2d and 3d [16]. 
The method presented here is independent of dimension 
and is also valid for fermions in the presence of exter- 
nal potentials. For a non-interacting fermion system the 
ground state is given by filling low energy modes ^i. The 
entropy of a region A may be expressed as [17] 

Se{L) = -tr[MlogM + {1-M) log(l - M)] (7) 

where M is a matrix with elements Mij =< (f)i\PA\4'j >, 
where Pa is projection on the region A. 

The distribution of fermion mimber may be ex- 
tracted from the moment generating function xW — 
J2Pae'^" = det(l - M + Me*^) [17]. Such functions 
appear in the theory of Quantum optics, where one stud- 
ies the distribution of detected photons. More recently 
the analog of x was introduced in the theory of quantum 
transport in mesoscopic systems and is referred to as the 
full counting statistics of fermions [18]. 

The cummulants of fermion number (or, equivalently 
Sz) in region A are given by derivatives of logx at A = 
0. Taking derivatives give expressions of the form <C 
{SS.r >= {-ird^x\x=o = Qn,it^MK 

Now trM' can be extracted from the cummulants by 
inverting Q, i.e. trM' = Q,"^ < (68^)" »• The ma- 
trix is lower triangular (i.e. Q^^ = for i > n) 
and may be calculated to any desired order. For exam- 
ple, trM =< SSz », trM^ =< > _ < SS^ », 
trM3 =< 6S, > -| < SS^ > +1 « SS^, >. 

We may write the entropy as the series, convergent 
whenever the eigenvalues of M are away from and 1 

SE = -j:n=2 7^^^iM"--^+i^-Mr)= (8) 

= -En=27^iQi,n«iSSzr» 

+ELi c^Qkji-'^y « i^szy »). 

For any finite matrix M this series converges as each 

term is smaller then dim M, it is easy to check that in 
fact Se < log 2 dim M as should be. One must note how- 
ever that convergence may be slow, depending on the 
eigenvalues of M close to unity or zero. Fortunately, all 
the terms in the series are positive, so the lower bound 
improves by adding more terms. In fact, the first non- 
vanishing contribution, namely n = 2 in the series cor- 
responds to the variance, and turns out to scale in the 
same way as entropy in the free fermion case [14]. 

Remarkably, the coefficients (Q~^);,n appearing are 
universal - i.e. will fit any non-interacting fermion the- 
ory, independent of its eigenmodes. One possibility is to 
measure cold fermionic atoms in a trap; by measuring 
the histogram of particle number in a certain region of 
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the trap, one could estimate the entanglement entropy 
using Eq. (9) regardless of the particular trap details. 

The intimate relation between bipartite entanglement 
and fluctuations of a conserved quantity is not restricted 
to non-interacting systems, but extends also to the 
strongly disordered spin- 1/2 XXZ model: 

H = ^ Ji (Sf • S\+i + Sf ■ SVi+^ + \St ■ S-i+i) (9) 

i 

with — 1/2 < A < 1 and the Ji are positive and randomly 
distributed. This interesting class of interacting random 
1-d theories was recently discussed in Ref. [19], where it 
was shown that the bipartite entanglement entropy of a 
segment of length L with the rest of the chain is Se (L) = 
i In 2 log2 L. The ground state of the Hamiltonians in Eq. 
(9) consists of a frozen liquid of valence bonds, which is 
referred to as the random singlet phase [20-22]. 

The Hamiltonian (9) only commutes with •Sf^^^j = 
^ Sf (note, however, that its ground state has a full 

i 

rotational symmetry). Therefore Sa ~ "^f the op- 

erator of choice for estimating the entanglement between 
part A and the rest of the chain. In the random singlet 
phase there arc two types of singlets: (a) Nab singlets 
connecting between A and B, (b) Naa + Nbb singlets 
connecting sites in A to other sites in A, or sites in B to 
other sites in B. As explained in Ref. [19], each singlet 
contributes 1 to the entropy, and therefore: Se = Nab- 
In addition, each singlet contributes 1/4 to the variance 
of Sa- In this case indeed we have the relation: 

Se = ^{1\S\) = Nab- (10) 

Note that Eq. (9) still applies in this case where the 
full counting statistics of the z-direction spin becomes 
X(A) = = det((l - M)e-*^/2 _^ Me}^!"^), with 

M now being a diagonal matrix with entries 1/2 for each 
singlet of type (a), and otherwise. 

Eq. (10) raises the question whether there is an even 
more direct relationship between variance of conserved, 
quantized, quantities and entanglement in the context of 
spin chains. We considered this for resonating valance 
bond (RVB) states of six spins with varying weights for 
each singlet configuration, and found that although in 
most of the Hilbert space varS\ < Se, a small region 
exists where this inequality is violated. This is associated 
with the formation of strong ferromagnetic correlations. 
It is possible that by putting a few restrictions on a RVB 
state one can prove a general relation. Nevertheless, the 
variance of a conserved quantity of a quantized object 
can be thought of as an ad- hoc measure of entanglement. 

The interest in entanglement is wide spread due to the 
prospect of engineering and controlling entangled states 
in which two [or more] microscopic objects, although 
separated by a macroscopic distance, display quantum 



correlations. This gives even more iirgency to under- 
standing how to measure entanglement quantitatively. 
To determine whether they have such exotic states, ex- 
perimenters carry out Bell inequality violation tests (c.f. 
Refs. [23, 24]). Could the schemes presented here real- 
istically quantify entanglement in many body systems? 
So far we considered the entanglement entropy in the 
ground state, thus implicitly assuming that the temper- 
ature is T = 0. Realistically, however, the temperature 
is finite, the state under consideration is not in a pure 
state, and Eq. (1) should, in principle, be generalized 
to account for thermal fluctuations. However, one can 
always imagine obtaining a lower bound on the ground 
state entanglement entropy by extrapolating the appro- 
priate measurement entropy, as detailed above, to T = 0. 
We leave the investigation of the efficacy of such realistic 
procedures for future work. 

We thank J. Preskill and L. B. loffe for discussions. 
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